The Bogoliubov Theory of a BEC in Particle Representation 
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In the number- conserving Bogoliubov theory of BEC the Bogoliubov transformation between 
quasiparticles and particles is nonlinear. We invert this nonlinear transformation and give general 
expression for eigenstates of the Bogoliubov Hamiltonian in particle representation. The particle 
representation unveils structure of a condensate multiparticle wavefunction. We give several exam- 
ples to illustrate the general formalism. 
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I. INTRODUCTION 

Experimental realization of a Bose-Einstein condensate 
(BEC) in trapped alkali atoms triggered feverish activ- 
ity in cold atom physics Q . Recent development of ex- 
perimental techniques in trapping and cooling of dilute 
atomic gases allows for preparation, manipulation and 
control of a condensate in various states and various trap- 
ping potentials. For instance, excitation of a condensate 
to vortex and solitonic states became possible |2(- Cold 
atoms are also becoming a convenient toy system to study 
solid state physics. For example, cold atoms trapped in 
an optical lattice potential undergo quantum phase tran- 
sition to Mott insulator phase when the barrier between 
minima of the lattice potential is increased 0, • 

Theoretical description of a many body system is a 
challenge to physicists. At low temperature our under- 
standing of BEC is usually based on a mean field approx- 
imation within the Gross-Pitaievskii equation (GPE) 
GPE provides satisfactory experimental predictions when 
depletion of a condensate, induced by interactions be- 
tween atoms, is negligible. In order to describe these 
quantum fluctuations around the condensate in a sys- 
tematic way one can use the Bogoliubov theory of BEC. 
The key idea of the original Bogoliubov theory [f| is the 
U(l) symmetry breaking approach where the atomic field 
operator is assumed to have a nonzero expectation value. 
This coherent state necessarily involves superposition of 
states with different numbers of atoms, an assumption 
very far from experimental reality. 

A number-conserving Bogoliubov theory has been in- 
troduced in 0, 0- While taking into account a finite 
number of atoms in a condensate, this theory represents 
eigenstates of the Bogoliubov Hamiltonian as Fock states 
with well defined numbers of Bogoliubov quasiparticles 
in all Bogoliubov modes. In this quasiparticle representa- 
tion it is possible, and even pleasant, to calculate a single 
particle density matrix. However, higher order correla- 
tors are increasingly more painfull to get. For all practi- 
cal purposes in the quasiparticle representation the mul- 
tiparticle correlations remain hidden. The single particle 
density function carries very limited information about a 
many-body system, as is best illustrated by an example 
considered by Javanainen and Yoo Q . They simulated a 



measurement of atom density distribution in an atomic 
Fock state \N,N) (it is a model counterpart of an ex- 
periment with interference of two separate BECs [Tol|). 
With two counter-propagating plane waves as the two 
modes in the Fock state the single particle density func- 
tion gives uniform density distribution. However, sam- 
pling of the many-body probability density shows that 
each realization of the experiment reveals an interference 
pattern. Only after averaging over many realizations one 
gets the uniform density distribution predicted from the 
single particle density function. 

It seems thus crucial to understand the structure of 
a condensate multiparticle wavefunction especially when 
the number of non-condensed atoms becomes consider- 
able and may play a significant role in an outcome of 
a measurement. In this paper we derive an expres- 
sion for Bogoliubov eigenstates in particle representation 
i.e. as superpositions over Fock states with well defined 
numbers of particles (atoms) in a condensate and non- 
condensate single particle states. For example, the state 
without any quasiparticles (Bogoliubov vacuum) has a 
pair-correlated form 



44 + Xr 



it at 



10) 



(1) 



m>0 



where the index m runs over a basis of modes (j> m (x) 
orthogonal to the condensate wavefunction 4>o(x). This 
pair-correlated ansatz has been proposed in a review by 
Leggett but (except for a homogeneous condensate 
[Tl|) the eigenvalues A m and the eigenmodes (j> m remain 
unknown. In this paper we construct a simple relation- 
ship between these eigenmodes/eigenvalues and the ele- 
mentary Bogoliubov excitations of the system for a gen- 
eral inhomogeneous condensate which is of current ex- 
perimental interest. The pair correlated vacuum state is 
a foundation on which one can build the theory of BEC 
entirely in terms of A^-particle wavef unctions. 

The paper is organized as follows. In Section [H] 
we briefly summarize the number-conserving Bogoliubov 
theory introduced by Castin and Dum |8( for a system 
with a finite number of particles. In Section HTTl we derive 
the particle representation of the Bogoliubov eigenstates. 
In Sections II VI VI Vll we illustrate the particle representa- 
tion with a series of examples: a double well, a triple well, 
and a condensate in a harmonic potential. We conclude 
in Section EIII 
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II. BOGOLIUBOV THEORY WITH A FINITE 
NUMBER OF PARTICLES 

The original Bogoliubov theory (BT) 6| was build 
around the concept of spontaneous symmetry breaking: 
the field annihilation operator ip acquires a large nonzero 
expectation value, (ip(x)) = (po(x). This nonzero mean 
is sometimes regarded as a symptom of Bose-Einstein 
condensation. The original BT is build by expansion 
in small quantum fluctuations around this large mean, 
ip(x) = 4>o(x) + 5ip(x). This method has to be used with 
a lot of care. For example, careful application of the 
BT leads to the divergence of the quantum fluctuations 
that is interpreted as a quantum phase spreading of the 
condensate [l2j. In the original spontaneous symmetry 
breaking approach the condensate state is not a station- 
ary state because for {ip{x)) ^ the number of particles 
is not well defined and states with different numbers of 
atoms have different energies. A self-consistent version 
of the Bogoliubov theory, which explicitly takes into ac- 
count that there is a well defined number of atoms in a 
trap, was put forward in Ref. Q and independently by 
Castin and Dum ||. In the present section we briefly 
summarize the latter results. 

In the formalism of Ref. [8j it is assumed that most 
atoms are condensed, i.e. they occupy the same conden- 
sate wavefunction <po{x). To the leading order the state 
of the system is a Fock state 



1^0 : N) , 



(2) 



with all N atoms in the condensate wavefunction <p$. In 
this Fock state 



b : N\ V l^o : N) = , 



(3) 



as expected for any state with a well defined number of 
atoms. Fluctuations around such a perfect condensate 
are calculated with the help of an expansion 



In dimensionless units a Hamiltonian for N interacting 
atoms in a trapping potential V(x) is 



H 



dx 



-V^V^ + V(x)ftip + -ftflipip 



(6) 



where g is proportional to the s-wave scattering length. 
The Hamiltonian © is expanded in powers of Sip. 
To zero order in dtp the <fio satisfies a stationary GPE 



- g v '0o + V(x)<f>o + Ng \<j> \ 2 <po = M0o • (7) 

For a (f>o that solves the GPE the first order term in 
the expansion of the Hamiltonian © in powers of 5tp is 
zero. For small depletion (small number of atoms de- 
pleted from the condensate wavefunction) and N ^S> 1 
the second order term can be approximately expressed in 
terms of a number-conserving operator 



As an example of this approximation take a term 

2g5^ (pQal4> a 6ip = 
|2 8$ %ao6 , ip 



(8) 



2gN\cPo\ 



N 



2gN\4> \ 2 tfA-2gN\cp 
2 ff A^ | 2 AtA . 



2 8ip f Stp 
N 



(9) 



The last approximate equality is justified for small deple- 
tion. 

After all such approximations the second order term 
in the expansion of the Hamiltonian becomes 



tp(x) = a <fi (x) + Sip(x) . 



(4) 



Here ao annihilates atoms in the condensate wavefunc- 
tion <po while the fluctuation operator Sip annihilates in 
all modes orthogonal to the condensate. To make sure 
that the <po in Eq. J3J is indeed the condensate wavefunc- 
tion it is further assumed that there is no coherence be- 
tween the condensed and non-condensed parts, 



at 8tp 







(5) 



This constraint ensures that (po is an eigenstate of the 
one-body density operator. 

Perturbative expansion. — The theory is build 
by expansion in powers of the "small" Sip. This opera- 
tor is "small" when density of depleted (non-condensed) 
atoms is small as compared to density of the condensate, 
(5^(x)5iP{x)) < N\<p (x)\ 2 - 



H 2 w 1 J d 3 x (A^ - A) C 



A 

At 



(10) 



where 



C 



and 



H, 



GP 



gNQ\(p \ 2 Q 



gNQcplQ* 
H GP -gNQ\M 2 Q 



-gNQ*<Pl 2 Q 



H GP = --V 2 + V(x) + Ng \fo\ 2 



(11) 



(12) 



Here Q is a projection operator on the subspace orthog- 
onal to (po , 



Q 



(13) 
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Bogoliubov modes. — In order to diagonalize C 
one can first solve the Bogoliubov-de Gennes equations 



- - V 2 C/ m + V{x)U m + 2g\ct) \ 2 U m + g<t>lV m = 

fJ-U m + LO m U m , 

-\v 2 V m + V(x)V m + 2g\<f> \ 2 V m + g{^ ) 2 Ur, 



The eigenmodes of C are 



(14) 



Vrr, 



(15) 



This projection is necessary because the operator C 
in Ea. i|ll|) is different from the differential operator in 
Ea. (|14f> . H2 is diagonalized by a Bogoliubov transforma- 
tion 



H%) = hrn u ™(x) + v* n (x) . (16) 

m 

Given the normalization conditions 

(U m \u m >) - (V m \v m >) = 5 m m' , (17) 

6 m 's satisfy bosonic commutation relations 



[b m , b m ,} = 5 mm > and [b m , b m -] = 



(18) 



These commutation relations are valid in the present or- 
der of the expansion i.e. for small depletion. The w m (b m ) 
operators create (annihilate) quasiparticles. The diago- 
nalized H2 is a sum over harmonic oscillators 



(19) 



Equations 1)18119(1 define the BT in the quasiparticle rep- 
resentation. 

The Bogoliubov vacuum is an eigenstate of Hb without 
any quasiparticles 



b m I Oj,) = , for all m 



(20) 



In Section lnil we derive the particle representation of |0&) 
state and other Bogoliubov eigenstates (i.e. the states 
that are a result of action of the creation operators b^ 
on the Bogoliubov vacuum) . 



III. BOGOLIUBOV EIGENSTATES IN 
PARTICLE REPRESENTATION 

For small depletion the approximate Hamiltonian is a 
sum of harmonic oscillators, see Ea. (|19[l . The eigenstate 
of the Hamiltonian l|19f) without any quasiparticles is the 
Bogoliubov vacuum, compare Eq. l|2U[) . The formal so- 
lution for I Ob) in the quasiparticle representation is not 
suitable for any analysis of the structure of the eigenstate 
like e.g. comparison with exact diagonalization for model 
systems. In the following we give an expression for the 
Bogoliubov eigenstates in the particle representation. It 
allows one to investigate how atoms are depleted from 
the condensate wavefunction when strength of interac- 
tion between particles increases. 

&'s in terms of a's. — Using the expansion in Ea. i|16|) 
and the orthogonality relations in Eq. (|17|l we can express 
&'s in terms of A, 



b m = (u m \A) - (v m |A t ) 



(21) 



With definitions 



.*,t 



{Um\H) 



(22) 



and Eq. JSJ) we get the nonlinear Bogoliubov transforma- 
tion 



ATi/2 



(4 



a m 



(23) 



ao annihilates particles in the condensate, while u rn and 
v m annihilate in the modes u m {x) and v m {x) which by 
construction are orthogonal to the condensate wavefunc- 
tion <po(x). 

We note that the right hand side of the Bogoliubov 
transformation between quasiparticles and particles 
is nonlinear in particle operators. This nonlinearity is 
necessary to conserve the total number of atoms. We 
will effectively invert this nonlinear transformation when 
the total number of particles N is even. 

The finite N Bogoliubov vacuum |0b : N). — To 
write I Ob) in particle representation we need a particle 
operator a 1 that commutes with all quasiparticle annihi- 
lation operators b m , 



bm. , d) 



, for all m . 



(24) 



It turns out that a two particle operator of the form 

d = a Q a Q + ^2 z ki Ofca/ (25) 
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is the smallest particle number solution of Eas. (|24(l . The 
double summation in the definition of d in Eq. (|25|l runs 
over an orthonormal basis {4>i(x), fo^x), . . .} in the sub- 
space orthogonal to the condensate wavefunction 4>o(x). 
With the vanishing commutator in Ea. (|24|l and for an 
even number of atoms N we can write down Bogoliubov 
vacuum in the form 



\0 b :N) oc (rft)f | ) 



(26) 



which is annihilated by every 6 m , 6 m |0f, : N) = 0. This 
particle representation of |0f,) has a finite and well defined 
number of particles N. 

Excited states of the Hamiltonian 119|) are created from 
| Of, : N) by application of the quasiparticle creation op- 
erators 



bt 



ATl/2 



(27) 



Equation (|26[) gives particle representation of the Bo- 
goliubov vacuum for even N. For an odd N there is no 
iV-particle state annihilated by & m 's in Eq. (|23|l . This is 
easy to check for N = 1. The ground state of the Bo- 
goliubov Hamiltonian is not annihilated by 6 m 's. This 
is a direct consequence of the fact that 6 m 's are only 
approximate bosonic operators: the bosonic commuta- 
tion relations l|18fl are fulfilled only approximately when 
depletion is small. For large iV there is no significant 
difference between N and N + 1 so one may apply the 
results of an even particle system to an odd one. In the 
following we will illustrate this problem with an exactly 
solvable example. 

Solution for the matrix Z. — To get an explicit 
| Of, : N) we must solve the condition lt2^|) with respect to 



= 



1 



b m , d f 



(28) 



N 1 / 2 V 
At 



A/-1/2 



«m) ,44 + E z kia k a 

kd^O 

2 «m + E Zkl (["™' 5 fe] 



a i + [u m ,ai\a k 



k,l^0 



It is convenient to define matrices U and V , 



When u m 's make a complete basis in the subspace or- 
thogonal to 4>o , then U m k is an invertible matrix, and we 
get 



z = u- 1 V 



(32) 



The choice of the actual basis {4>i, 4>2, ■ ■ .} is in principle 
arbitrary. 

Diagonalization of Z. — We note that the matrix 
Z introduced in Ea. (|25|l is symmetric 

Z = Z T . (33) 
Under an additional assumption that Z is real, 



Z 



Z* 



(34) 



the matrix Z can be diagonalized in an orthonormal 
eigenbasis {<f>i, 4>2, ■ ■ ■} with real eigenvalues A^, and an- 
nihilation operators {d\, &2, ■ ■ .}, 



E 

k^Q 



\k a k 



(35) 



The diagonal d gives the finite- N Bogoliubov vacuum 
in particle representation in the form, see Ea. H26|) . 



|0 6 : N) 

N/2 

E 

no,Tii,."=0 



(Jt)f| ) 



(36) 



V(2n )!...(2n M )! ^ 

°N,2n + ...+2nAr : : x 



n ! 



,n M \ 



A? 1 ...\™\2no,...,2n M ) • 



Here for the sake of convenience we truncate to M Bo- 
goliubov modes. The ket \2n , 2ni, . . . , 2um) is a Fock 
state with 2no, 2ni, . . . , 2um particles in the eigenbasis 
{</>o, <j>i, ■ ■ . , 4>m }■ When necessary the particle represen- 
tation (|36|l can be translated into an iV-particle wave- 
function. 

As we can see in Ea. H36f) . every eigenmode can be oc- 
cupied only by an even number of particles. Particles are 
depleted from the condensate in pairs (2ni + . . . + 2um 
is even) , and every non-condensate eigenmode <pk can be 
occupied only by an even number of atoms. 

Single particle density matrix. — The eigenstates 
{</>o, (/>i,...} of the matrix Z turn out to be also eigen- 
states of the single particle density matrix, 



[u m , a k ] 



{Um\<Pk) 



vjn = V mk a k 



(29) 
(30) 



Using the completeness of the basis {</>i, </>2, ■ ■ ■} wc 
rewrite the last line of Eq. 1)28(1 as a matrix equation 



V = u z 



(31) 



pV>(x,y) = (0 b :N\ft(x)$(y)\Q b :N) 



OO 

E 

k=0 



(Of, : N\a< k a k \0 b : N) ^)My) 



(37) 



Off-diagonal elements vanish thanks to the even occu- 
pation numbers in the state l|36|l . Ea. (|37|l provides a 
familiar interpretation of the eigenmodes {4>q,4>\, . . .}. 
The next Section opens our series of examples. 
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FIG. 1: Number of depleted atoms from a condensate state 
as a function of frequency Q [see Eq. 1381 1. The solid line 
corresponds to the exact results, the dashed line to Eq. 149 . 
while the dotted line is depletion calculated directly from the 
Bogoliubov eigenstate in the particle representation, Eq. I|5i-Sfl . 
The inset shows the overlap of the exact state and the Bo- 
goliubov state |0{, : JV), Eq. 15211 . as a function of Q. Data for 
the double well system with N = 100 atoms. Q ~ 0.01 is the 
crossover point between the Mott insulator and superfiuid. 
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FIG. 2: Ground state \tp) of a condensate in a double well 
potential for N = 100 atoms and for Q = 10 (a), Q = 0.001 
(b). Circles (connected with a dashed line) correspond to 
the exact results while crosses are related to the Bogoliubov 
eigenstates in the particle representation. The ground states 
\ip) are projected on Fock states \N — n, n) where N — n is a 
number of atoms in the left well and n is a number of atoms 
in the right well. 



IV. 



DOUBLE WELL 



Let us put N atoms in a double well potential. In the 
tight binding approximation, when we restrict only to the 
Hilbert space spanned by the ground states in each well, 
this system is described by a boson Hubbard model 



ci + c 2 



V2 



(41) 



The only mode orthogonal to 4>q is 



2\V 



= - n 



c\c 2 



J'=l,2 



Here 



(38) 

for j — 1,2 is an operator of the number 
of atoms in the j-th well. The hjirij — 1) terms describe 
repulsive interactions between atoms in each well. The 
first term in the Hamiltonian is responsible for tunneling 
between different wells. All units were rescaled so that 
the only dimensionless parameter is the tunneling rate O. 

On the level of the GPE this system is described by a 
condensate wavefunction which is a vector of two complex 
amplitudes, 



with an annihilation operator 



oi 



c\ - c 2 



V2 



(42) 



(43) 



The double well system has been studied theoretically 
|13| and its multiwell generalizations were a subject of 
recent experiments Q. The ground state properties of 
this system are well known. In the regime of £1 ^> N the 
ground state is a quasi-coherent state 



( 



'0,1 



50,2 



(39) 



one amplitude for each well. Solution of the GPE for the 
ground state results in a symmetric state 



00 = ^(1,1). 
Atoms in this condensate are annihilated by 



(40) 



(4+4)^10) ~ (« ) JV |0) , (44) 

where all N atoms are in the same condensate wavefunc- 
tion (|40|l and the number of depleted atoms N\ — (a\a\) 
is zero. 

In the Mott insulator regime of Nil <C 1 the ground 
state is a Mott insulator, 



N 



|0) . 



(45) 
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This state is very far from the condensate JUJl, it has a 
huge fraction of depleted atoms ^ = ^. One half of 



all atoms are depleted from the condensate wavefunction 
(HI predicted by the GPE. 

Quasiparticle representation. — The Hubbard 
model (|38ll has been studied in the framework of the BT 
only very recently [lj] . The Bogoliubov-de Gennes equa- 
tions are 



N 
N 

— Q(m — U 2 ) + — (u 2 + V 2 ) = +UJU2 , 

N 

-Q(V2 - Vl) + y («i + Ui) = -UJVi , 

N 

—Cl(vi - v 2 ) + y (^2 + u 2 ) = -uov 2 , (46) 

The index j = 1,2 in Uj and Vj numbers a well. These 
equations give only one Bogoliubov mode orthogonal to 
the condensate <po 



X 



-1 



VA 2 - 1 



(47) 



where 
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FIG. 3: The first excited state \tp) of a condensate in a double 
well potential for iV = 100 atoms and for Q = 10 (a), f2 = 
0.001 (b). Circles (connected with a dashed line) correspond 
to the exact results while crosses are related to the Bogoliubov 
eigenstate with one quasiparticle 6' [0b : N). The states \ip) 
are projected on | TV — n,n) states where N — n is a number 
of atoms in the left well and n is a number of atoms in the 
right one. 



X 




P. 2 



(48) 



In the quasiparticle representation the number of de- 
pleted atoms in the Bogoliubov vacuum state is approx- 
imately given by 



dx AU = 



dx S^clq— alStp 



{V V 



X 2 -l 



(49) 



The approximate equality is justified when the depletion 
is small and N ^> 1. This dN is plotted with a dashed 
line in Fig^ When il — > 0, then the dN diverges like 
dN ~ fi -1 / 2 . As N is finite, the divergence in dN is 
unphysical. The divergence of dN is not present when 
we calculate depletion directly from the Bogoliubov state 
in the particle representation. 

Particle representation. — The operator b for the 
Bogoliubov mode (|T7I) is 



HA) - HAt) 



XciqCii + d\do 
V^(A 2 - 1) 



(50) 



We look for a two-atom annihilation operator d ~ dg + 
Aidf such that b commutes with d\ [b,tv] = 0. The 
solution is 



X 



(51) 



This operator defines the Bogoliubov vacuum in the par- 
ticle representation 



|0 6 : N) ~ (<?) 2 |0) . (52) 
The number of depleted atoms in this state, 



dN= (0 b : AT|atai|0 b : N) 



(53) 



is plotted with a dotted line in FigQ] It compares surpris- 
ingly well with dN in the exact ground state (solid line) 
even in the Mott insulator regime (Nfl « 1). In con- 
trast, in the insulator regime the standard Bogoliubov 
theory (dashed line) gives an unphysically divergent dN. 

In Fig. [21 we compare the Bogoliubov eigenstates 
with the exact ones in both weak and strong tunneling 
regimes. We find surprising agreement between both so- 
lutions that is present even in the Mott insulator regime. 
To see how the Bogoliubov vacuum smoothly interpolates 
between the quasi-coherent state l|44|) for large £1 (large 
X), and the Mott insulator gSJ) for small Q (X — ► 1) it is 
convenient to rewrite Eq. I|51|l as a product of annihilation 
operators in two in general non-orthogonal modes, 



fx 



ao + a - . (54) 



Ol 

X 



7 



With this d the Bogoliubov vacuum becomes a quasi- Fock 
state 



|0„ : N) 




(55) 



The two modes become the same for large fi (large X), 
and they become orthogonal for Cl — ► (X — ► 1). This 
quasi- Fock representation may be useful in any case when 
one can approximately truncate the Hilbert space to only 
one Bogoliubov mode. For the double well this trunca- 
tion is exact. 

In Fig. |21 we compare the first excited state of the dou- 
ble well system with the Bogoliubov eigenstate with one 
quasiparticle b'\Qb : TV) for two values of il, one in the 
Mott insulator and the other in the superfluid regime. 
Like for the ground state we find surprisingly good agree- 
ment between the two states even in the Mott regime. 

The surprising agreement between the Bogoliubov 
eigenstates and the exact eigenstates that we found even 
in the Mott insulator regime is rather an exception than 
a rule, as we show for a triple well system. 



V. PERIODIC TRIPLE WELL 

The Bogoliubov theory applied to the double well po- 
tential turns out to be surprisingly good. In the present 
example it does not work so well and we can investi- 
gate how the Bogoliubov approximation stops working 
with increasing depletion from the condensate. The pe- 
riodic triple well system is described by a boson Hubbard 
Hamiltonian 



H 



:S\Y 



1 3 

n E + oE • (56) 



We proceed similarly as for the double well. The ground 
state of the GPE is 



6o = -^ (1,1,1) 



(57) 



and it is associated with an annihilation operator &Q. 
There are two modes orthogonal to <T>o 



1 
1 

71 



, 27Ti 27T1 

l,e + 3 , e 3 



l,e 3 ,e + 3 



(58) 



with annihilation operators a + and d_ respectively. The 
Bogoliubov-de Gennes equations give two Bogoliubov 
modes, let us call them + and — , with 
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FIG. 4: Ground state of a condensate in a triple well po- 
tential for TV = 12 atoms and for Q = 4 (a), Q = 0.1 (b). 
Circles correspond to the exact results while crosses to the 
Bogoliubov eigenstates in the particle representation. The 
states are projected on |no,n+,n_) states where no is a 
number of atoms in a condensate wavefunction Eq. I57H while 
n+ and n- correspond to numbers of atoms in the orthogonal 
modes c/>+ and <A_ Eq. 1581 . If the depletion is small (in (a) 
dN ~ 0.03) the Bogoliubov eigenstates match well the exact 
ones. However, for a large depletion (in (b) dN « 1.9 that 
is significant as compared to TV = 12) the exact states be- 
come considerably different from the Bogoliubov prediction 
that only states \N — 2i,i,i) should give nonzero contribu- 
tions. Panel (c) shows the overlap of the exact and Bogoli- 
ubov eigenstates as a function of Q for 57 £ [0.02, 4]. 
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FIG. 5: Exact ground state of a condensate in a triple well 
potential for TV = 13 atoms and for Q = 4 (a), Q = 0.1 (b). 
The states are projected on |no,n+,n_) states where no 
is a number of atoms in a condensate wavefunction Eq. I57H 
while n+ and n_ correspond to numbers of atoms in the or- 
thogonal modes <f>+ and </>_ Eq. 1581 . Contributions of states 
that correspond to even numbers of depleted atoms from the 
condensate, i.e. |TV — 2i, ...... .), are marked by full circles. 

Compare the structure of the eigenstates with the similar re- 
sults for TV = 12 presented in Fig. 01 
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X<f>± -<f>± 
u+ = = , v± = . = 

Vx 2 -i vx 2 -i 



(59) 



where 



X= 1 




These modes give 



b± = 



Xa[ ) a± + aoati 



VN(X* - 1) 
The operators b± commute with a cF such that 



d = a a — — a+a_ , 
which can be easily diagonalized by operators 



(61) 



(62) 
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FIG. 6: Eigenvalues Xj of Z matrix corresponding to the 
ground state of a condensate in a ID harmonic trap for 
gN = 20 (open circles) and gN = 50 (full circles). Note 
that with increasing gN more and more modes get significant 
contributions to a condensate eigenstate, see Eq. 1361 . 



Oi 



a + + a_ 



72 



«2 = 



V2 



with eigenvalues 



Ai — — A2 — — 



1 

X 



(63) 



(64) 



In the Bogoliubov theory only pairs of atoms are de- 
pleted from a condensate wavefunction. Recently, in ex- 
act stochastic calculations Carusotto and Castin also 
observed such a pairwise depletion from a condensate in a 
harmonic trap. In Fig.0^ we show the structure of a Bo- 
goliubov vacuum together with the exact ground state of 
the triple well system which confirms the predicted pair- 
wise depletion. However, with decreasing tunneling rate 
fi, that makes the interaction between atoms relatively 
stronger, the exact ground state looses the pairwise struc- 
ture expected in the linear Bogoliubov theory, see Fig.^J). 
This is a clear signature of a breakdown of the linearized 
Bogoliubov theory. 

In Fig. 03 we show an exact ground state for an odd 
number of particles. Again like for even N particles are 
depleted in pairs. Although we were not able to find 
any simple form for the Bogoliubov ground state with 
odd TV, this example suggest that, especially for large N, 
accurate predictions can be extrapolated from even N. 
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FIG. 7: Plots of a condensate mode 4>o(x) and three the most 
significant eigenmodes <f>j{x) to which the condensate is de- 
pleted for the system in a ground state of a ID harmonic trap 
for gN = 50. Note that (f>j(x)'s are eigenmodes of p fl '(x,y), 
compare Eq. l37ll . 



Bogoliubov vacuum for ground and solitonic states of a 
condensate in a ID harmonic trap 



VI. BOSE-EINSTEIN CONDENSATE IN A 
HARMONIC TRAP 

In the previous two exactly solvable examples the con- 
struction of the finite N Bogoliubov vacuum was rela- 
tively easy because there were only one or two Bogoli- 
ubov modes. In this Section we analyze the finite N 



V(x) = -x 2 . (65) 

The solitonic state is the first excited antisymmetric state 
oftheGPE0j. 

We constructed the Bogoliubov vacuum in particle rep- 
resentation in the following steps. To begin with we 
found the relevant stationary state 0o of the GPE in 
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FIG. 8: Eigenvalues of Z matrix corresponding to a soli- 
tonic state of a condensate in a ID harmonic trap for gN = 50. 
Compare the present plot with the results for a condensate in 
a ground state presented in Fig. [(J 
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FIG. 9: Plots of a condensate mode (j>o(x) and three the most 
significant modes <j)j( x ) to which the condensate is depleted 
for the system in a solitonic state of a ID harmonic trap for 
gN = 50. Note that ^ 3 -(a;)'s are eigenmodes of p' 1 '(x,y), 
compare Eg . 1371 . 



a harmonic potential: the ground state or the solitonic 
state. Then we diagonalized the Bogoliubov-de Gennes 
equations (|14|) for a given background </>o and projected 
the obtained modes on the subspace orthogonal to the 
condensate wavefunction like in Ea. (|15(l . In this way we 
obtained the Uj and Vj modes of the number-conserving 
Bogoliubov theory. In the next step we constructed ma- 
trices U and V like in Eas. (|29l3U[) . In order to construct 
such matrices we had to choose a basis in the subspace 
orthogonal to the condensate wavefunction cf>Q. To make 
the matrices finite we had to truncate the number of Bo- 
goliubov modes to a finite M. We got fast convergence 



of results with increasing M for a basis constructed out 
of the uj modes, 

4> 2 ~ u 2 - 0i(0i|«a) , •■■ ■ (66) 

Finally, we solved Eq. (|31[1 with respect to Z and then di- 
agonalized the obtained Z matrix. The eigenvalues were 
ordered with decreasing modulus, |Ai| > |A 2 | > ... > 
| Am |- As is clear from the state Infill , even a small de- 
crease of | Afc | from one mode to the next results in a sub- 
stantial drop in the average number of atoms occupying 
this mode. 

In Fig. we plot the spectrum of the Z matrix, cor- 
responding to the ground state of the condensate in a 
harmonic trap for two different interaction strengths gN. 
With increasing gN more and more modes have signifi- 
cant contribution to the condensate eigenstate. Because 
pairs of atoms are depleted from the condensate, the sin- 
gle particle eigenmodes 4>k are allowed to be both even 
and odd as shown in Fig.0 Note that <j>j are also eigen- 
modes of the single particle density matrix, see Eq.J27|). 

For a condensate in a solitonic state, in the spectrum of 
the Z matrix (Fig. |HJ the dominant eigenvalue Ai is pos- 
itive (contrary to other eigenvalues) and the correspond- 
ing anomalous mode is localized in the soliton notch |l6| , 
see Fig. [5] The anomalous mode <f>i (x) is a (nearly) zero 
mode related to translational motion of the dark soliton 
with respect to the condensate. 



VII. CONCLUSION 

In the present paper we have constructed the Bogoli- 
ubov eigenstates in the particle representation. The par- 
ticle representation gives much more intuitive insight into 
physical processes that are responsible for depletion of 
a condensate. Comparison of the derived eigenstates 
with exact solutions for model systems reveals that as far 
as a number of depleted atoms is small the Bogoliubov 
eigenstates match well exact eigenstates. In the Bogoli- 
ubov theory the condensate is depleted by pairs of atoms. 
However, when strength of the interaction between atoms 
increases the linear Bogoliubov theory breaks down and 
the condensate starts being depleted by odd numbers of 
particles. 

We have applied our approach to analyze structure of 
a realistic condensate in the ground and solitonic states 
of a harmonic trap. Having calculated eigenstates in the 
particle representation one can in principle extract any 
information about a system and simulate results of any 
measurement. Examples of such realistic applications 
will be given in our future publication |17| . 
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